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ABSTRACT 

We investigate the stability of partially ionised, differentially rotating, diffusive disc 
threaded by both azimuthal and vertical magnetic field. The general stability crite- 
rion of such a disc in the presence of axisymmetric fluctuations can be stated purely 
in terms of ambipolar and Hall diffusivities. It is shown that the disc is magnetorota- 
tionally unstable if the sum of scaled ambipolar and Hall diffusivities are larger than 
some numerical constant determined by the rotation profile of the disc. This criterion 
■ suggests that subthermal diffusive discs are always unstable to almost radial fluctu- 

ations. The field geometry and obliqueness of wavevector (encapsulated together in 
the topological factor g), plays dual role of not only assisting MRI in ambipolar dom- 
inated disc but also making otherwise stable region in Hall-ambipolar diffusion plane 
unstable. 

The vertical magnetic field and transverse Alfven fluctuations are fundamental to 
magnetorotaional instability in Hall-Ohm diffusion dominated discs. The disc having 
both azimuthal and vertical field in the presence of oblique axisymmetric fluctuations 
can be mapped to a purely vertical fleld threaded disc with purely vertical wavenum- 
ber. Although, azimuthal field and obliqueness of wavevector rescales the growth rate 
and window of instability, such a scaling permits MRI to operate much closer to the 
midplane in Hall-Ohm dominated disc than is otherwise possible for a purely vertical 
field threaded disc. 

The maximum MRI growth rate in the ambipolar diffusion dominated disc is 
proportional to the topological factor g. For long wavelength fluctuations the maximum 
growth rate is two fifth of ideal MHD case (3/4 The short wavelength fluctuations 
are likely to be damped by the diffusion. 

The diffusive discs are prone to new kind of MRI-the diffusive MRI which has no 
counterpart in non-diffusive discs. The diffusive MRI is caused by interplay between 
advection and diffusion of the field. However, since differential rotation of the disc 
is also at the heart of diffusive MRI, the maximum growth rate of the instability is 
identical to ideal MHD case. 

The excessively diffusive accretion discs in addition can also become unstable due 
to interplay between differential rotation and field diffusion. The ensuing diffusive 
instability grows at same maximum rate as MRI. However, the transport of angular 
momentum in the disc may not be efficient due to diffusion instability. 

Key words: accretion discs, plasma, Stars:Formation, MRI 



1 INTRODUCTION 

It was Fred Hoyle (1960) who first suggested that only 
a magnetohydrodynamic model can explain the outward 
transport of angular momentum and resulting slow rota- 
tion of the protosolar nebula. Although, Hoyle's proposal 
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was well received by the scientific community (Safronov 
1972), the detailed mechanism of magnetic field-assisted 
angular momentum transport remained elusive for another 
three decades before Balbus & Hawley (1991) discovered 
such a mechanism: the magnetorotational instability (MRI). 
Though the stability of magnetised Coutte flow was analysed 
by Velikhov (1959) and Chandrasekhar (1961), their work 
remained largely unnoticed until Balbus & Hawley (1991) 
not only clarified the underlying physical mechanism of this 
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instability but by specifically applying it to accretion discs, 
showed that this instability will play an important role in 
angular momentum transport in the disc. Only after Bal- 
bus & Hawley's work did it become clear that a weak, sub- 
thermal magnetic field can facilitate the local transport of 
angular momentum in a weakly magnetised Keplerian disc. 
The most attractive feature of this instability is that the 
growth rate depends only on the differential rotation of the 
disc, an abundant source of free energy. The identification 
of this powerful local instability has provided a firm phys- 
ical foundation for the paramctcrisation of viscosity in the 
turbulent discs (Shakura & Sunyaev 1973; Hawley & Balbus 
1991; Stone et al. 1996; Balbus & Hawley 1998; Balbus 2003, 
2009). 

The initial formulation of MRI by Balbus & Hawley 

(1991) was for a fully ionised, ideal inagnetohydrodynamic 
(MHD) medium in which field and fluid are well coupled 
to each other, a far cry from the solar nebula where often 
nonideal MHD effects such as Ohm, Hall and ambipolar dif- 
fusion, play an important role. In fact, a large fraction of 
matter in the universe is not in a fully ionised state but 
is partially ionised with varying degree of ionisation deter- 
mined mainly by local thermodynamic conditions. It is well 
known that in partially ionised discs viz. the solar nebula, 
protoplanetary discs, the discs around cataclysmic variables 
and black holes, the number of plasma particles, particularly 
away from thermal source, is often very small compared to 
the neutrals and thus, the plasma particles are not frozen 
in the magnetic fleld but can slip through it due to col- 
lisions with neutrals (Mestel & Spizer 1956). This realisa- 
tion led to the generalisation of MRI to weakly ionised discs 
(Blaes & Balbus 1994) where collision between the neutral 
and plasma particles are major source of energy dissipation. 
However, if the neutral-ion collision time is faster than the 
rotation period of the disc, the medium moves as a single 
fluid and MRI in such a medium remains unaffected by col- 
lisions. In the opposite limit, when the neutral-ion collision 
time becomes comparable or exceeds the rotation period, 
resistive quenching of MRI occurs in the disc (MacLow et 
al. 1995; Jin 1996; Hawley & Stone 1998; Sano & Miyama 
1999; Turner et al. 2007; Bai & Stone 2011). 

The magnetic fleld in a differentially rotating diffusive 
disc can evolve either (a) due to advection by the fluid (when 
diffusion is weak), or, (b) due to both advection and diffu- 
sion, or, (c) only due to diffusion (when advection by the 
fluid is negligible). Whereas advection dominated discs are 
described in the ideal MHD framework, remaining cases, (b) 
and (c) are important when non-ideal MHD effects can not 
be neglected. When (b) is valid, discs can become MRI un- 
stable along one of the two possible pathways: (i) magnetic 
fluctuations are time dependent and its evolution is both 
due to diffusion and fluid advection (Wardle 1999; Balbus 
& Terquem 2001; Kunz & Balbus 2004; Desch 2004; Wardle 
& Salmeron 2011) and, (ii) magnetic fluctuations are almost 
time independent and fleld drift is exactly offset by the ad- 
vection. We shall name this last case, diffusive MRI (DMRI) 
(Wardle &: Salmeron 2011). However, when differentially ro- 
tating disc is excessively diffusive that is when diffusion dom- 
inates advection of the fleld, there is no MRI or, DMRI but a 
new kind of diffusive instability (DI) appears in the disc even 
when there is no fluid motion to transfer the angular mo- 
mentum (Rudiger& Kitchatinov 2005; Kunz 2008; Bejarano 



et al. 2011; Wardle & Salmeron 2011). Therefore, presence 
of favourable magnetic field topology in diffusive discs opens 
up new channels through which shear energy can be trans- 
ferred to waves. 

The realisation that non-ideal MHD effects not only 
causes dissipation but can also assist MRI emerged only af- 
ter it was shown that Hall effect, depending on the sign 
of magnetic liclicity , can significantly modify the instabil- 
ity (Wardle, 1999, hereafter W99, Balbus & Terquem, 2001, 
hereafter BTOl). The role of ambipolar diffusion in assist- 
ing MRI was soon discovered (Kunz & Balbus 2004; Desch 
2004). It was shown that the ambipolar diffusion (AD) not 
only causes anisotropic dissipation, but under favourable 
field geometry can also assist the instability. More impor- 
tantly, unlike ideal MHD where angular velocity must de- 
crease outward in the disc. Hall and AD modified MRI is 
not dependent on the angular velocity profile (Balbus & 
Terquem 2001; Kunz & Balbus 2004; Desch 2004). Only the 
magnetic field topology and direction of wavevector becomes 
central to the MRI in the partially ionised discs (Kunz & 
Balbus 2004; Desch 2004) and role of angular velocity pro- 
file is relegated to insignificance. All in all, non-ideal MHD 
effects plays dual role in the disc: not only it causes dissipa- 
tion, but can as well assist MRI in the presence of favourable 
magnetic field topology. 

Often, astrophysical discs are not uniformly ionised. 
The radial as well as vertical profile of the ionised mat- 
ter depends on the proximity of the disc to the thermal 
and non-thermal sources. For example, the surface layer of 
the protoplanetary discs (PPDs) is ionised by both cosmic 
(Umcbayashi & Nakano 1981), and x-rays (Glassgold et al. 
1997; Igea & Glassgold 1999; Feigelson 2005; Ercolano et al. 
2008). However, close to midplane, matter is largely shielded 
from external ionising sources and thus, except for the tiny 
ionisation by radionuclides (Umcbayashi & Nakano 2009), 
matter is mostly neutral. As a result vertical structure of the 
disc is layered, with surface layer well coupled to the mag- 
netic field, whereas midplane of the disc, magnetically dead 
(Gammie 1996; Stone et al. 2000; Balbus 2003; Salmeron 
& Wardle 2003, 2008; Terquem 2008; Bai & Goodman 2009; 
Turner et al. 2007, 2010; Armitage 2011). Clearly, the nature 
of magnetic diffusion through the disc changes with chang- 
ing ionisation profile. In PPDs, Ohm-Hall is the main diffu- 
sion closer to the midplane whereas ambipolar-Hall diffusion 
dominates the disc surface (Wardle 2007; Bai 2011). Since 
MRI in such a disc can not operate closer to the midplane, it 
raises an interesting question: can local magnetic field topol- 
ogy at least extend the domain in which MRI operates, par- 
ticularly closer to the disc midplane? As has been shown 
recently by Wardle & Salmeron (2011; hereafter WSll) the 
answer to this question is in the affirmative. They showed 
that if the vertical field is aligned to the rotation axis, the 
MRI can be excited at much lower ionisation level than is 
otherwise possible. Therefore, ambient ionisation together 
with the local field topology determines the efficiency of MRI 
in transporting angular momentum particularly close to the 
midplane of PPDs. 

In the present work we investigate the MRI in a dif- 
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ferentially rotating, partially ionised disc threaded by both 
azimuthal and vertical magnetic fields, i.e. B = (0, B^, Bz). 
Wc shall assume that the fluctuations are axisymmetric and 
the wavevector k has radial and vertical components, i.e. 
fe = {kr,Q,kz)- It is pertinent to note here that the stabil- 
ity of such a diffusive disc was analysed by Kunz & Balbus 
(2004) and Desch (2004). Notwithstanding the simultaneous 
appearance of these works, the general dispersion relation 
was analysed only by Desch (2004) who not only provides a 
clear explanation of the underlying physical mechanism of 
MRI in ambipolar diffusion dominated discs but also identi- 
fies crucial topological parameter g (which encapsulates disc 
field topology and obUqueness of wavevector) that assists 
the instability. However, as we shall see, analytical results 
of Desch (2004) pertains to excessively diffusive discs, that 
is when fluid advection of the fleld is unimportant. Clearly, 
several important questions remain unanswered. For exam- 
ple, is it possible to state a general stability criterion for 
diffusive discs? What is the relative role of fluid advection 
and field diffusion in destabilising the disc? Does ambipolar 
diffusion affects the window over which Hall diffusion causes 
MRI? Can we identify the most MRI favoured field geometry 
in the disc ;' How does the maximum growth rate in various 
diffusive regimes compares with the ideal MHD case? How 
is the critical wavelength of the instability affected by the 
diffusion and the field? To clarify these issues, and to com- 
pare the behaviour of a magnetorotationally unstable disc in 
various diffusive regimes, we reinvestigate the stability of a 
diffusive disc in the present work. Thus, notwithstanding the 
similarity of our formulation with above mentioned works, 
not only the treatment of the problem but also the set of 
questions that we wish to address here are quite different. 

The paper is organised in the following fashion. In sec- 
tion 2 we briefly describe the basic set of equations suitable 
for the magnetised, non self gravitating, diffusive disc ro- 
tating around a central point mass. In section 2.1 we give 
linearised equations in terms of diffusion tensor before de- 
scribing general dispersion relation. By using a proper scal- 
ing, the dispersion relation is reduced to a very simple form 
which allows us to infer the general stability criterion in 
diffusion dominated discs. In following subsections, we show 
that diffusion plane, depending on most unstable and cut-off 
wavenumbers can be subdivided in three distinct regions. In 
section 3, dispersion relation is analysed in two dimensional 
Hall-Ohm (subsection 3.1) and Hall-ambipolar (subsection 
3.2) diffusion planes. In section 4, limiting cases of dispersion 
relation is discussed. Section 5 contains several subsections 
discussing possible applications of results to PPDs. A brief 
summary of results is presented is section 6. 



2 FORMULATION 

We shall start with standard description of a non self- 
gravitating, differentially rotating, magnetised, thin disc of 
partially ionised matter which is threaded by both toroidal 
and poloidal magnetic fleld B — (0, B^, Bz). We adopt cylin- 
drical (r , , z) coordinate system anchored at the central 
star of mass M and assume that the angular velocity of 
the disc, f2(r) is constant on the cylinder. The velocity of 
the fluid in the local frame corotating with the disc angu- 
lar velocity f2(r) can be written as v = V — Vti>{ro) 4>- Here 



^0(^0) = rQ,{rQ), V is the velocity in the laboratory frame 
and f2(ro) is fiducial angular velocity at ro. The time deriva- 
tive in the local frame dt becomes dt + ^(ro) in the labo- 
ratory frame. The continuity and momentum equations are 



|+V.(p.) = 0, 



(1) 



^ + {v ■V)v-2nv^r+^vr4> = - — + ^^^ , (2) 
ot 2 Si p cp 



where the epicyclic frequency x is 



2 2 dfi^ 
X =40 -I- 



d\nR 



(3) 



Various terms in the momentum equation have their usual 
meaning. The induction equation is 



dB 
'dt 



Vx 



{v X B) — J ■ 



4-rrrjH 



J xb 



+^ (J X 6) X 6 



(4) 



where b is the unit vector along the magnetic field. The 
expression for Hall, r?H, ambipolar, tja and. Ohm, r)o 
for weakly-ionised ion-electron-neutral plasma are (Konigl 
1989; Balbus & Terquem 2001) 



cB 
4n eue 



DB^ 



Vo 



where 



7j 



C me'JePn 

47r , 



rrii + rUn 



(5) 



(6) 



and < a V >j is the rate coefflciont of collisional momentum 
exchange. Above diffusivities can be written in compact form 
as (Pandey & Wardle 2008) 

rjH= (^) ,7]A = D (^) ,and,r,o=P:'vH, (7) 
\ojH J yvni y 

where D = p„/p is the ratio of the neutral mass density 
to the bulk density, p = pn + Pi, va = B/y^A-rr p is the 
Alfven velocity, /3e = tOce/ven is the electron Hall parameter, 
which is a ratio of the electron cyclotron to electron-neutral 
collision frequencies. The Hall frequency is ujh = piCOd/p- 
Although, above diffusivity expressions are for a dustless 
medium, our analysis is quite general and is equally appli- 
cable to partially ionised medium containing dust. One could 
simply use a more general expressions for diffusivities, e.g. 
Wardle (2007). 

Current J is given by Ampere's law. 



J= -^VxS. 

47r 



(8) 



The role of magnetic field diffusion in assisting or in- 
hibiting MRI can be made transparent if, following WSll, 
we recast the induction Eq. (4) as 



dB 

dt 



= Vx 



(^ + ^,)xB-i^j|i 



(9) 
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where the magnetic field diffusion or, drift velocity'^ 



vb = riP 



(V X B)^ X b 



- Vh 



(V X B), 



B 



(10) 



and rjp — rjo + rjA is the Pedersen diffusivity. The parallel 
and perpendicular components of the current in the preced- 
ing equations refer to the orientation with respect to the 
background magnetic field. Note that first term on the right 
hand side of Eq. (9), in the bracket V Xv X B is the usual 
advection term. In the absence of diffusion, when vb ~ 
and rjo = field is advected by the fluid at finite velocity 
V. The second term Vxvs X B in the bracket is due to 
slippage or, drift of the field linos through the fluid at finite 
velocity vg. In the fiuid frame, the magnetic drift velocity 
Vb contains the effect of field diffusion. 

In the absence of parallel current, i.e. setting last term 
to zero in Eq. (9), the induction equation reduces to famil- 
iar form except now unlike ideal MHD, flux ^ — f B ■ dS 
through the arbitrary surface S encircled by a contour C in 
the plasma is not frozen in the fluid frame but in a frame 
moving with v + vb- Such a formulation of induction equar- 
tion allows us to visualise magnetic field as real physical en- 
tity that is drifting through the fluid with a given velocity. 
In this frozen-in frame, with the field and fiuid moving at 
V + Vb, we have cE± + {v + vb) X B = 0, which in the 
presence of (VxB)|| becomes 

f = ,/^(VxB)„.d5. (11) 

It is clear from Eq. (11) that the presence of parallel current 
and not Ohm diffusion is responsible for the flux nonconser- 
vation in a weakly ionised plasma. Since energy dissipation 
is inevitable in such collisional plasmas, at flrst sight this 
result appears paradoxical. The non-ideal MHD effects, par- 
ticularly Pedersen (Ohm -I- Ambipolar) diffusion, owing to 
its dissipative nature, must affect magnetic flux. However, 
energy loss comes from the redistribution or, relaxation of 
magnetic field and not from its annihilation (e.g. magnetic 
reconnection) (Parker 1963). Therefore, total flux is con- 
served even in the presence of Ohm diffusion only if parallel 
current is absent in the medium. 



2.1 Dispersion relation 

We shall assume a thin differentially rotating disc with 
V ,Cs ,va <S vk which is threaded by a magnetic field 
having both azimuthal and vertical components, i.e. B = 
{0, B^, Bz). Further, a general rotation profile f2(r) oc r"'' 
is assumed. The radial variation in the physical quanti- 
ties is neglected by assuming krr S> 1. Although non- 
axisymmetric fluctuations can play an important role in 
magnetic fleld amplification in the disc (Balbus & Hawley 
1992a; Terquem & Papaloizu 1996), our analysis is limited 
to axisymmetric fluctuations only. Thus, linearising Eqs. (1) 
- (4), and assuming axisymmetric perturbations of the form 
exp {ik ■ X + at), with k = {kr ,0 ,kz), the continuity equa- 
tion (in the Boussinesq approximation) becomes 



k-6v = 0. 



(12) 



^ We shall use diffusion and drift terms interchangeably as they 
both mean same physical process. 



The (r , 0) components of the momentum equation be- 
comes 



Sv ■■ 



i k fi 



2ki 
a 



[a^ + 2kl (2-g)] V (? - 2) 
and magnetic field drift velocity can be written as 



SB, (13) 



5VB = — 



6z r?p sf)ii-\- g ffp 



5B, 



(14) 



where 5v = Sv/va, Svp ~ 5vb/va, SB — SB/B, k — 
kvA/i^, o = a/Q., b = B/B, kz = kz/k, n = ■ bj = kzbz 
and f) = 7) 0./v\ . The topological factors g and helicity s are 

g = —kr kzb^bz , s = filcz = (k ■ b) (k-Cl) . (15) 



Here f2 = n/O. In table 1 we provide the list of important 
symbols and their usual meaning to ease reading. The im- 
portant role of geometric factor g in assisting AD modified 
MRI has been noted by Desch (2004). We shall see that not 
only maximum growth rate but also parameter window of 
the AD assisted MRI depends on this geometric factor. 

The linearised induction equation (r , 4> components) 
becomes 



+ 



[a2 + 2fei (2-g)] 



+k' 



(5-2) 

r]r<t, 



2k^ 
a 



dB = 0, (16) 



where the diffusivity tensor have following components 
Vrr = rjo +blrjA , r}r4, = s'nn+ gfiA, 
V<l>r = {gVA- stIh) /kl , rj4,4, = f)o + - fcr flA ■ (17) 
The dispersion relation obtained from Eqs. (12)- (16) 



a* + P {rjp + frr) + C^a^ + Cia + Co = . 
where 

C2 = fe"* (r)p TjT + fji) - qk^ (s fjH+ gVA) 
+2 ki (2 - g) + 2 P , 
Ci = P [2pz {2-q) + fc'] iVP + Vt) , 
;^ P -I- ((4 - q) S17H - q gVA) P -2q fc^j k" 

-2kzP q (2 - g) {sr]ii+ gVA) 
+2 P kl (2 - q) {fjp fjr + m) 



(18) 



Co 



and 



(19) 



(20) 



The dispersion relation, Eq. (18) is well known in the liter- 
ature (Kunz & Balbus 2004; Desch 2004). However, as has 
been noted above, only Desch (2004) analysed the general 
dispersion relation. Wc shall re examine above dispersion 
relation de novo to answer following questions. What is the 
general criterion under which a diffusive disc is magnetoro- 
tationally unstable? Does ambipolar diffusion affect the do- 
main over which Hall destabilises the disc? What is the role 
of g in Hall diffusion driven MRI? How important is the 
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Table 1. List of frequently used symbols. 



Sybmol 


Explanation 


Symbol 


Explanation 


V 


advection velocity 


Vb 


r lelcl driit velocity 


V 


Velocity 


k 




hz 


kz/k 


bz 


Bz/B 


hj. 
^(p 


j^(p / ^ 






s = ^kz 


Helicity 


g = -kr kz bz 


Topological factor 


g' 


gjkz 


b 


B/B 


c 




cr 


a/kz 


k 


kVA/^ 


k 


kbz 


VP =VO+ VA 


Pedersen Diffusivity 


Vt = VO + VA 


Total diffusivity 


flp (vt) 


VP (vt) 


fip (Vt) 


»)p Wt)/ (p.bz) 


ifjo,flA) 


(vo,riA) 


(fio,fiA) 


ivo ,flA) / it^bz) 


m 


7]H ^/v\ 


VH 


m/bz 



role of vertical field and vertical wavevector in exciting MRI 
in diffusive discs? How does the presence of azimuthal field 
affect the instability? Is unstable wavevector window a func- 
tion of disc magnetic field? Docs favourable field topology 
in a diffusive disc also implies that MRI can operate much 
closer to the midplane than the limit set by purely vertical 
field threaded disc (e.g. WSll)? To reiterate, our analysis 
has a different focus and is complementary to Desch (2004). 
Rescaling physical quantities as 



a— ^ a = a/kz, k k = kbz ,fip fjp = 



kzbl 



m m 



VP 

= 1^, (21) 
bz 



the dispersion relation, Eq. (18) can be recast in the follow- 
ing simple form 



ak + bk^ + c - 



where the coefficients a, b and c are 

c 



0, 



(22) 



a 



■ P — y , b = Sa ■ 



(23) 



with 



/3 = (»7h + VP Vt) 5-^ -|- {fjp + vt) a 
1 



+2 (2 - q) 



VPVfT -\- \VtA + 



S={i)p+rfr) [2 (2-g)+a'] + {^^] 



7 = 9 [9 VA+nti + 



(24) 



and g' = g/kz- The dispersion relation, Eq. (22) can also be 
recast as 



^¥ + 5k'^a+ [2 {2-q) + a'^] = 
7F [2 {2-q)+a^ 



(25) 



which shows that for q < 2, coefHcients S, /3 and 
[2 (2 — q) + a^^ on the left hand side are all positive. Thus 
positive a will require positive 7 as square bracket on the 
right hand side is always positive for q < 2. Therefore, a 



differentially rotating disc is unstable only if 



> 0. 



<1 [9 Va+'Vh + 



(26) 



We have arrived at the most general stability criterion for a 
diffusive disc: whenever (i) < g < 2 and (ii) g' va + Va + 
1/ (2 — g) > 0, a differentially rotating disc is unstable. Note 
that Ohm diffusivity does not appear in the above expression 
which is not surprising considering that Eq. (26) is the long 
wavelength limit of the sufficient stability condition, Co < 0. 

The physical condition for the onset of instability can 
be presented as a competition between the fluid advection 
and magnetic field drift in the plasma. This can be seen 
from the azimuthal and radial components of Eqs. (13) and 
(14) respectively, which suggest that when a = 0, the radial 
advection of the fluid, Svr is equal and opposite to the radial 
drift velocity 5vg^. This can be seen by combining azimuthal 
component of Eqs. (13) and radial component of Eq. (14)] 



5vr + Svbt = ik h 



VA SBr 



(27) 



Here we have neglected Ohm diffusion. By setting 7 = 
near marginal state, in the vicinity of t;^ <C 1, we get Svr + 
5vg^ Rs 0. This provides a simple physical explanation for 
why diffusion will assist MRI. The outward radial slippage of 
the field due to diffusion weakens magnetic tension force. As 
a result magnetic restoring force in the wave gets dominated 
by the inertia force which causes increased inward radial 
drift of the fluid element. This is how diffusion assists MRI. 

As has been noted in the introduction, the magnetic 
diffusion opens up completely new channels through which 
shear energy can be transferred to waves. Therefore, it is nat- 
ural to ask whether above stability criterion, Eq. (26) could 
have been as well derived from reduced induction equation, 
i.e. dropping one or other term in Eq. (16). It will be shown 
later (please see section 3) that in the limiting case, when 
advection and diffusion terms in Eq. (16) balance each other, 
we arrive at the identical criterion. 

For a diffusive Keplerian disc, since q = 3/2, condition 
(i) q £ (0, 2) is automatically satisfied. Thus the stability 
criterion for a Keplerian disc becomes 



g'VA + m + 2>0, 



(28) 
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which states that the sum of ambipolar and Hall diffusiv- 
ities with suitable coefficients determined by the oblique- 
ness of wavevector and field topology should exceed some 
given numerical value predetermined by the Keplerian pro- 
file. Eq. (28) provides a particularly simple stability test for 
diffusive discs. For example, Hall-Ohm dominated disc is 
unstable if rjH > —2. This stability criterion can be stated 
in terms of Hall and rotation frequencies as 



f2 > — 2 6^ ujh = 2uiHZ- 



(29) 



where lohz is the Hall frequency defined in term of ver- 
tical field. For negative the disc is magnetorotationally 
unstable if Hall frequency is less than one half of the Keple- 
rian frequency. The Hall unstable region [e.g. Fig. 5, WSll 
or, Fig. 1-3, BTOl] in a purely vertical field threaded disc 
[bz = 1) becomes arbitrarily small in 6^ — >■ limit when the 
disc is threaded by non-vertical field. We see from Eq. (28) 
that at riH^/vA = — 2 6z, since tja is always positive, am- 
bipolar diffusion can assist MRI only if, (; > 0, in agreement 
with Desch (2004). 

From stability criterion Eq. (28) we see that whereas 
Hall-Ohm stable discs {r]H ^ — 2&zU^/n) may or may 
not be ambipolar unstable, Hall-Ohm unstable discs {riH > 
—2 6z v^/^) can always become ambipolar unstable for non- 
vertical field and oblique wavevector since g > can be eas- 
ily satisfied. The ambipolar diffusion, on the other hand, not 
only assists MRI when g > 0, but also enlarges parameter 
space over which Hall can assist MRI. Therefore Hall assists 
MRI when 



rjH > -2 6z va/0, 



■VA- 



(30) 



Since bz ,b^ £ [0, 1] and maximum g = 0.25, above inequal- 
ity implies that for non-zero rjA, Hall diffusion causes MRI 
at much larger negative value than for purely vertical field 
(WSll, BTOl). 

To summarise, stability criterion of a magnetised dif- 
fusive disc can be stated purely in terms of disc diffusivi- 
ties. The MRI in the Hall-Ohm dominated disc is assisted 
by Hall diffusion if its value exceeds some numerical value 
predetermined by the disc rotation profile. In ambipolar- 
Ohm dominated discs, in addition, the topological factor g 
(which carries the signature of non-vertical field and oblique 
wavevector) should be positive. The presence of such a pos- 
itive topological factor in turn enlarges the parameter space 
over which Hall assists MRI. The topological factor g plays 
the dual role: it not only assists MRI in ambipolar regime 
but also enlarges the domain of Hall assisted MRI. 

2.2 Division of diffusion plane 

It is clear from Eq. (22) that the wavenumber range over 
which MRI occurs extends between zero and the cut-off 
wavenumber 

1/2 



2 (2 - g) 7 


2 (2 - g) 


ijpfn + (vK + 


-7 



(31) 



The disc is magnetorotationally stable if 7 0. In the op- 
posite limit, when the disc is unstable, the diffusion space, 
like WSll can be divided into three different regions de- 
pending on the denominator in Eq. (31). When (a) < 7 < 
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^ III 








k„-> - 
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Figure 1. The division of diffusion plane according to tlie values 
of wavenumbers fee , fco . Region I in diffusion plane correspond to 
finite fcc , fco- Region II correspond to finite kc finite and fcc — ^ 00 
and region III correspond to {kc , fco) oo- 



2 (2 



1 

2-q 



the cut-off wavelength is 



finite [Fig. (1) labelled I]. On the other hand, when (b) 

. 2' 

f/p fjT + [fjn + ■ 



1 

2-q 



the cut-off wavenum- 



7 ^ 2 (2 - g) 

ber occurs at infinity. The unstable region in the diffusion 
plane in this case can be further subdivided depending on 
the value of most unstable wavenumber kg — —bjla which 
is 



Kq 



[2 (2-g) + ag] 7 - ^ arp 
2/3-7 



(32) 



Here a"o is maximum growth rate which can be derived by 
setting 6^ - 4 a c = in Eq. (22). 

Thus unstable region when fee — > cjo is divided in region 
II and III depending on whether 7 > 5 (Tq/ [2 (2 — a) -f ctq] 
or, 7 < (55o/ [2 (2 - a) -f- jg] [Fig. (1) labelled II and HI]. 
Region II correspond to finite fco and infinite kc whereas re- 
gion HI correspond to infinite fco and kc. We shall utilise this 
analysis in the following section to examine the dispersion 
relation. 



3 ANALYSIS OF DISPERSION RELATION IN 
TWO DIMENSIONAL DIFFUSION PLANE 

Coupling between the magnetic field and matter varies with 
height in protoplanetary discs. The detailed calculation of 
diffusivities at 1 AU for a minimum mass solar nebula in 
which dust has settled to the midplane of the disc suggests 
that if the disc is threaded by a pG field, then Ohm is the 
dominant diffusion from midplane to three scale heights^ 
at 1 AU whereas Hall dominates between three and five 
scale heights (Wardle 2007; Bai 2011). The ambipolar dif- 
fusion is dominant beyond five scale heights. The presence 
of dust does not change this picture either though at 5 AU, 
Hall becomes dominant diffusion very close to the midplane 
[Figs. (10) and (11) in Wardle (2007) or Fig. 2 in Bai (2011)]. 
However, for stronger magnetic field \B £ (0.5, 2) Gauss], 

which is determined by the hydrostatic balance between gravity 
and pressure gradient 
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when dust has settled to the midplane, at lAU Hall is 
the dominant diffusion mechanism between midplane and 
5 scale heights. In the presence of 0.1 fim dust, not Hall but 
Ohm becomes dominant between midplane and two scale 
heights (Wardle 2007; Bai 2011). We may conclude that 
with changing scale height, as one moves from midplane 
to surface layer in PPDs, diffusivity typically changes from 
Ohm to ambipolar with Hall sandwiched in between, i.e. 
{rjo riH ^ tia) with a large overlap region where Hall- 
Ohm or, Hall-ambipolar is simultaneously important. There- 
fore, owing to the layered structure of PPDs it is impor- 
tant to carryout the analysis in two dimensional, Hall-Ohm 
iVH — rjo) or, Hall-ambipolar {r]H — 'riA) planes. Furthermore, 
such an analysis provides a clear physical interpretation of 
the dispersion relation. To keep the analysis transparent we 
restrict ourselves to a Keplerian disc and assume q = 3/2. 



a 0.5 

o 
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Stable 


1 

X II Untable 







-3 -2 -1 -0.5 0.5 
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-3 -2 -1 -0.21 0.5 



3.1 Hall-Ohm diffusion plane 

It is clear from above discussion that owing to its layered 
structure. Ohm and Hall is dominant diffusion close to the 
midplane of the protoplanetary disc. Therefore, we shall first 
analyse the dispersion relation, Eq. (22) by retaining Ohm 
and Hall diffusion and set rjA = 0. The reduced dispersion 
relation in two dimensional diffusion planes provides a sim- 
ple pictorial representation. This can be seen by writing the 
dispersion relation, Eq. (22) as 



1 + 



where 



4 



A ■- 



1 1 

1+^2 + J2 



3_A 

4 



(33) 



(34) 



Eq. (33) is identical to Eq. (BIO) of WSU except now the 
physical quantities are normalised according to Eq. (21). Re- 
call that in WSll, disc is threaded only by vertical magnetic 
field and wavevector is parallel to the magnetic field. That 
the dispersion relation derived for a disc with non-vertical 
field and waves propagating at an oblique angle can be 
mapped to a purely vertical field with vertical wavenumber 
suggests that the vertical field and transverse fluctuations 
are fundamental to MRI in the Hall-Ohm diffusion domi- 
nated discs. In the absence of either vertical field or, verti- 
cal wavenumber disc is stable to axisymmetric fluctuations. 
Therefore, we may conclude that the Hall-Ohm dominated 
disc is magnetorotationally stable if waves are purely radial 
or, ambient field is purely azimuthal. We stress that this 
conclusion is valid only for the axisymmetric fluctuations. 
However, disc threaded by a purely toroidal fleld do become 
magnetorotationally unstable in the presence of nonaxisym- 
metric fluctuations (Balbus & Hawley 1992a; Terquem & 
Papaloizu 1996). 

The above dispersion relation, Eq. (33) describes a cir- 
cular locus in riH — rjo plane for given o and k. Fig. (2) 
shows the dependence of the marginal locus on bz- Recall 
that magnetorotationally unstable region in a vertical fleld 
threaded disc with k = (0,0, fc^) correspond to ryu > —2 
[Eq. (28); see also Fig. 5 of WSll or, Fig. 1-3 of BTOl where 
vs. 2r)H has been plotted]. Since in this case we see from 
Eq. (27) that fm > -2 (6^ = L = 1) implies that Hall 
assists MRI only when the radial drift of the field exceeds 



Figure 2. The solution of Eq. (33) in the Hall-Ohm diffusive 
plane is shown in the above figure. In Fig. 2(a), purely verti- 
cal magnetic field {B^ = B) and fc^ = 1 is assumed whereas in 
Fig. 2(b), presence of both azimuthal and vertical fields (Bz = 
0.5 i?) is assumed keeping unchanged. The vertical line di- 
vides stable region from magnetorotationally unstable region. The 
Hall assisted MRI region in the presence of purely vertical field 
[6z = 1 , J7H > -2, Fig. 2(a)] shifts to the right [fju > -2 6^ = -1, 
Fig. 2(b)] for non-vertical field. 



radial advection of the field, i.e. Sv^^ > Svr- Thus, the in- 
stability interval (—2,-1/2) in which all wavenumbers are 
unstable correspond to 5v^^ > Svr and Sv^^ < Sv^f, since 
form Eqs. (13) and (14) we get 



5v,p = ——ik SBr 



and, Sv^ 



-i k rjH 5Br 



(35) 



When j)h > —1/2, both radial and azimuthal drift of the 
field exceeds corresponding drift due to fluid advection. In 
this case, large range of wavenumbers in region I (finite fcc 
and fco) in Fig. (2) are destabilised. 

To summarise, when only radial field drift 5v^^ exceeds 
radial advection of the fluid 5vr, magnetorotationally unsta- 
ble waves of arbitrary kc and fco (corresponding to all three 
region I, II and HI) in Fig. (2) is assisted by Hall diffusion. 
However, when both radial and azimuthal fleld drift exceed 
respective component of the fluid advection velocity, unsta- 
ble waves with only flnite kc and fco (corresponding to region 
I) is assisted by Hall diffusion. 

The maximum growth rate in the Hall-Ohm diffusion 
plane can be easily inferred by noting that the radius of 
the circle is determined by the non-zero shifted part of the 
?7o-axis in Fig (2). Thus setting 770 = in Eq. (33) and 
equating the shifted part to the right hand side of Eq. (33), 
i.e. {BAf = {2,A/Af yield 

ao = ^|fc-n|. (36) 

Note that above maximum MRI growth rate in diffusive 
Hall-Ohm dominated disc is identical to the fully ionised, 
ideal MHD discs. This is not surprising since behaviour 
of the fluid in Hall-MHD is similar to the ideal MHD ex- 
cept that instead of magnetic flux freezing, generalized flux 
J (ujh b -I- V X d) • dS is frozen in the fluid (Pandey & Wardle 
2006, 2008). We note that maximum growth rate depends 
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Figure 3. The solution of Eq. (38) for kz = 1 and ao = 0.5 Q. 
Contours of ag = 0.5 r2 (straight Unes) for fez = 1 (sohd) and bz = 
0.5 (dotted) is shown in the figure. The corresponding semicircles 
Eq. (33), are the boundaries between regions I and II (like in 
Fig. (2)) for fe^ = 1 (solid) and = 0.5 (dotted) and fc oo is 
also shown in the figure. 



The most unstable wavenumber becomes 



(2 (To rjo - 3/2 sryn) + 



M^(-3fci+2ag) 



(39) 



When vertical field is weak and Ohm diffusion is negligible, 
ko can become arbitrarily small. However, limit on the field 
is set by vertical scale height H ~ Cs/fl of the disc owing to 
the requirement fcoH > 1. 

The parameter window of the instability, is between in- 
terval [0, kc] where Eq. (31) for Hall-Ohm case becomes 



2bz 



Vo + ('?H + 2bz 



(i?H + I 6z 



(40) 

Region I encompasses both ideal-MHD as well as dif- 
fusive limits. The inner boundary of region H occurs when 
kg ^ oo and ctq is given by Eq. (38). Following expression 
for the locus separating regions H and HI is derived from 
WSll [Eq. (B15)] 



2klhz i^kl + ASrl 



{kl+ai) ( 



16 



(41) 



on the orientation of the wavevector which is also similar to 
the ideal MHD case (Balbus & Hawley 1992b). 

As has been noted above, the rjH — rjo plane can be di- 
vided into three distinct regions which traces three distinct 
form of (T as a function of k^ . In region I, both most unsta- 
ble wavenumber, ko and cut-off wavenumber, kc are finite 
whereas in region II corresponding to fyn G [—2,-1/2], fco 
is finite and fcc — ^ oo. In region III, both fco and kc. — )■ cxd. 
Rescaling Eq. (Bll) of WSll, region I, in Fig. (2) outside 
the semi-elliptical locus is described by 



-2 , 2 

vo+ ^^ 



9 2,2 

16 



(37) 



with eccentricity y^l — ^^//i. 

The maximum growth rate derived by rescaling 
Eq. (B13) of WSll yield 



24 o-o 



9fc2 bz 



16 b: 



■ Vo 



2kz bz 
kl + &l 



(38) 



and as shown in Fig. 6 by WSll, the contours of constant 
6"o are straight lines in region I. The orientation of the field 
and wavevector merely modifies the slope of the curve. 

It is clear from Eq. (36) that MRI growth rate is maxi- 
mum when kz = 1. We can also infer from Eq. (38) that only 
vertical field threaded disc (bz = 1) grows at maximum rate. 
In order to see this, we plot Eq. (38) for constant ao with 
kz = 1 m r)o — rjH plane for bz = 1 (solid line) and bz = 0.5 
(dashed line). Making use of Eq. (33), we also show corre- 
sponding contours for above values of field which subdivides 
the diffusion plane between finite and infinite k. We see that 
when bz = I, maximum growth rate line has a higher slope 
in riH — rjo plane than when bz — 0.5. This implies that MRI 
will grow at maximum rate only when the disc is threaded 
by purely vertical field. Therefore, MRI grows at maximum 
rate in discs having only vertical field and waves propagating 
along the field (kz — 1). 



rjo = 



kl bl ao (dkl 



16 a^ 



(fci-Fag) (9fci + 16a2) 



(42) 



which traces an arc in fju — fjo plane from fjo — fju ~ 
{0,—2bz) to fjo — fjn = (0, — |6z) as ao runs from to 
0.75 fcz. Therefore, with decreasing vertical field, region HI 
becomes smaller (Fig. 2). 

We conclude from above analysis that Hall-Ohm diffu- 
sion dominated discs with very simple field geometry, having 
only vertical magnetic field and vertical wavevector plays 
fundamental role in destabilising it against axisymmetric 
fiuctuations. Since MRI grows at maximum rate in such a 
disc, we may infer that B = (0,0, Bz) and fc = (0,0, fcz) 
is the most magnetorotationally favoured topology in the 
disc. The role of the azimuthal field and radial wavevector 
is limited to scaling various regions of instability as shown 
in Fig. (2). 



3.2 Hall-Ambipolar diffusion plane 

The ambipolar diffusion is the dominant diffusive mecha- 
nism in the surface layer of the protoplanetary disc followed 
by the Hall in the middle layers. This generic picture is 
broadly valid close to the protosolar nebula as well as fur- 
ther out in the disc at 10 AU even when grains have not 
settled to the midplane of the disc (Wardle 2007; Bai 2011). 
Therefore, it is important to analyse the dispersion relation 
in the Hall-Ambipolar diffusion dominated disc. 

By setting 770 = the dispersion relation, Eq. (22) can 
be written as 



'2 I '2 i\r 

Vh + Va = N 



1 1 

l + a^^k^ 



(43) 
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where 
and 

1 



5/4 
1 + 5-2 



3/4 

fc2 



VA + 



2^2 



1 



(44) 



AT = 



2\2 



16 



^ ) 



(45) 



This dispersion relation describes a circular locus in rjA ~ fjH 
plane for given a and k. 

In order to find the maximum growth rate in ambipolar 
diffusion dominated disc, we set rjo ~ rjH ~ in Eq. (18) 
and treat it as an equation in u}\ = v\, and in 

ri plane. Taking the partial derivative in u>\, and k^ and 
eliminating oj^ from one of the equations, we arrive at the 
following equation for the maximum growth rate 

fjA I (l — /i^) P fjAd'o — (f + M^) (^5 + 3 3' fc^ 77A j CTQ 

-0.25 (i6m'- 9/) P,^a + y3'|=0. (46) 

Although a general expression for the maximum growth rate 
can be easily written from above equation, it acquires par- 
ticularly simple form in the limiting cases. In the long wave- 
length hmit (fc — s> 0) the maximum growth rate becomes 

3gn 



2 (l + M^ 



(47) 



It is not surprising that the maximum growth rate in am- 
bipolar diffusion dominated disc is proportional to g which 
encapsulates the topofogy of the disc field and obliqueness 
of the wavevector (Desch 2004). The parameter g couples 
the poloidal wavevector to the toroidal magnetic field. Since 
the rotational velocity = r Q is also along the azimuthal 
direction, the presence of nonzero azimuthal field along 
allows the fluctuations to tap into the reservoir of free en- 
ergy. Therefore, topological factor g acts like a switch to 
the reservoir and if this factor is positive, then anisotropic 
ambipolar diffusion is able to assist the instability. 

Recall that the stability condition, Eq. (26) which is 
long wavelength limit of the sufficient condition, Co < 
requires positive g. As g £ [—0.25,0.25], we conclude that 
the maximum growth rate in purely ambipolar regime is 
less than half of the ideal MHD or, Hall-MHD regime for 
H = 1/2. Since anisotropic dissipation is responsible for MRI 
in ambipolar diffusion dominated discs (Desch 2004), part 
of free energy, otherwise available to waves in nondissipative 
discs is lost to dissipation causing smaller MRI growth rate. 

The maximum growth rate in short wavelength (A: — >■ 
00) becomes 



3 (1 + M')flf^ , , 



3 (1 + m')5 



(48) 



where A = 27 g -I- 9 g -|- 16 fc^ (l - /x ) . Dropping terms 
of the order of {g'^ ■, to can be written as 



3 {\^ yi)gVL 



2s 



(49) 



2(1-M^) V(I-m') 
which shows that short wavelength fluctuations are largely 
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Figure 4. The solution of Eq. (43) for fcz = 0.3, 6z = 0.9 are 

shown in the top [g > (fe^ < 0)] and bottom [g < (fe^ > 0)] 
panels respectively. The vertical line separates stable and unsta- 
ble regions. The unstable region is subdivided in three different 
regions corresponding to different fco a-nd fcc as shown in Fig. (1). 



damped. For example when g — 0.25 and kz = = l/\/2, 
(JO = -0.19^2. 

Since growth rate in the ambipolar regime is propor- 
tional to g we may conclude that predominantly Alfvenic 
(/i 1) or, predominantly magnetosonic {fi — > 1) fluctu- 
ations will grow at very small rate since g — )■ 0. However, 
oblique waves with k^ = kr = l/\/2 and hz = — 6^ = l/v'2 
grows at maximum rate. Therefore, oblique magnetic field 
inclined at 7r/4 in 0z plane is most favoured MRI geometry 
in ambipolar diffusion dominated discs. 

In Fig. (4) the locus of circle for both positive and neg- 
ative g is shown. The vertical line separates the stable re- 
gion from the unstable region. The sign of g is changed by 
changing the sign of azimuthal magnetic field. It is clear 
from Fig. 4(a)-(b) that when <? > 0, area under the curve in 
TjH — rjA plane is larger than when g < 0. This happens be- 
cause ambipolar diffusion is conducive to MRI for positive g. 
Note that on the horizontal Hall axis, if ryn > —2bz— g rjA/s, 
the Hall diffusion will assist MRI. Therefore, for positive 
g and 77 a 7^ 0, Hall diffusion can affect MRI over regions 
which were otherwise inaccessible in a purely vertical field 
threaded disc (Fig. 5, WSll or, Fig. 1-3, BTOl). Therefore, 
when g > 0, ambipolar diffusion not only assists MRI but 
simultaneously allows Hall diffusion to operate in enlarged 
parameter space. As a result the locus dividing region I and 
II in Fig. 4(a) has an outward bulge. Similarly, negative g, 
?7a 7^ causes shrinkage of parameter space for Hall diffusion 
[Fig. 4(b)]. Clearly, favourable field topology and obliqueness 
of wavevector not only affects MRI in ambipolar diffusion 
dominated disc but also extend the domain over which Hall 
can operate in the disc. Therefore, the sign of g plays crucial 
role in affecting the stability of Hall-ambipolar dominated 
discs. 

In Fig. 5 we plot locus of the circle dividing regions I 
and II shown in Fig. (4) for various k^. We see that with 
decreasing k^, area under the curve increases. This implies 
that almost radial fluctuations of arbitrary wavelength is un- 
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Figure 5. The dependence of locus on fcz diving region I and II 
[ Fig. (4)] for 6z = 0.9 are sliown in the top [g > (b^ < 0)] and 
bottom [g < (6^ > 0)] panels respectively. 



stable in AD dominated disc, in conformity with the general 
criterion Eq. (28) or, Eq. (75). Although accessible riA — tjh 
parameter window increases with decreasing fcz, the growth 
rate of the instability, which is proportional to g oc fcz be- 
comes small. When g < 0, radius of the semi-circle locus 
decreases with decreasing kz, implying that almost radial 
fluctuations are severely damped. 

Like Hall-Ohm case, Hall-ambipolar plane can also be 
divided into three distinct regions. In region I, outside the 
semi-circular locus 



_ 9 
4.(1^^) ^ 16 



/ 2 , l2 



(50) 

describes a half-circle in the rjH ~ tja plane. Region I en- 
compasses both ideal MHD and ambipolar-Hall space. All 
the wavenumbers within region II inside the semicircle, de- 
scribed by Eq. (50) are unstable. 

The maximum growth rate in r/A — tjh plane can be 
found by setting the discriminant of Eq. (22) to zero. This 
gives 



6z Ri/k , 



Z A 

'7A 



Here 



Ri = 



hi 



S (qi 



ai I 



[11 



ai 



bi 



11.2 



(ai + ci) =F \/{ai -cif +4bl 



and 



ai =bl (9 - 16 (T,^) , 6i = 6/i (l + ctq) do , 



ci = 4 fc^ (1 + o-o)^ (d,i - 4 /i^ o-q) 



(51) 



(52) 



(53) 



(54) 



(55) 



with do = (l + /^^) (Jo — I g' ■ If B = Bz and wave is prop- 
agating along the field, then Eq. (51) reduces to Eq. (B13) 
of WSll since 



Ri 



24 ao (l + gg) 
(9-16a2) 



(56) 



The wavenumber corresponding to the maximum growth 
rate is —b/2a which is found by setting the discriminant 
b"^ - 4ac = in Eq. (22). 

In region I, both fco and kc are finite whereas in region II, 
fco is finite and fcc — >■ oo. Region HI in Fig. (4) is approached 
asymptotically when both most unstable wavelength, and 
critical wavelength approach infinity. The dispersion relation 
(43) in k va — >■ oo limit becomes 



77H + 



5/4 
1 + 5^ 



+ fl ??A + 



do 



2^2 (1 + *,^ 



iV 



2\2 



(57) 



The locus that separates region II and III, can be found 
by equating Eqs. (51) and (57). This gives the following 
parametric solution 



4fiRi ai/{i + a^) 



3/xiii-2A;2 (1 + 52) do i + 



4 Lao 



ibzRi~2kz (l + *fl) rfo 



(58) 



(59) 



which traces an arc from (—2 6^,0) to [rjH ,TiA)\ao) in the 
diffusive plane. 

We see from Eq. (31) that the wavenumber below the 
critical value 



kc VA 



T?A + (ryn + 2 6z 



('7H 



, (60) 



are unstable in the Hall-ambipolar disc. Only when the 
wavevector is aligned to the ambient field, i.e. = 1, Eq. (60) 
becomes identical to Eq. (40) if we replace r/A with rjo ■ How- 
ever, even for a purely vertical field, Eq. (60) does not re- 
duce to Eq. (40). Note that when wavevector is aligned to 
the ambient field, i.e. = 1, Ohm and ambipolar diffu- 
sion can be combined together as Pedersen diffusion and 
their effect on the wave is identical-they both cause damp- 
ing (WSll). However, when /i 7^ 1, whereas Ohm diffusion 
isotropically damps the waves, damping by ambipolar diffu- 
sion is anisotropic. Therefore, /i 7^ 1, lifts the diffusive de- 
generacy causing ambipolar diffusion to assist MRI. Thus, 
critical wavelength in this case can not be reduced to Hall- 
Ohm case even for a purely vertical field. 

To conclude, anisotropic nature of ambipolar diffusion 
which is encapsulated in the topological factor g, plays cru- 
cial role in assisting MRI in Hall-ambipolar diffusion dom- 
inated discs. Not only positive g assists MRI in ambipolar 
regime but simultaneously enlarges the parameter space over 
which Hall can influence MRI. In a purely ambipolar diffu- 
sion dominated disc, maximum growth rate is proportional 
to ig/4 when the wave is almost vertical whereas it is 3 (?/2 
for almost radial waves. 
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4 LIMITING CASES OF THE DISPERSION 
RELATION 

In the standard MRI picture pertaining to fully ionised, dif- 
ferentially rotating disc radial field is generated by the in- 
ward motion of the fluid element and the rotational shear 
generates azimuthal magnetic field (SB^) from radial field 
(SBr) (Balbus & Hawley 1991). In a partially ionised, dif- 
fusive disc, not only differential rotation but both Hall and 
anisotropic ambipolar diffusion are also capable of generat- 
ing azimuthal field from the radial field. Therefore, magnetic 
diffusion provides additional new pathways through which 
free shear energy can be transferred to waves. Even when one 
of these pathways (through which free energy is channelled 
to waves) is blocked, disc can still become magnetorotation- 
ally unstable in the presence of diffusion. To see this, we first 
note that three terms in Eq. (16) are ~ 1 rj and thus, 

one of the following three scenarios may exist in a diffusive 
disc (WSll). 

A. Ideal MHD:'^ Diffusion in the disc is weak in com- 
parison with the fluid advection, i.e. field is frozen m the 
fluid; ^5v ~ dsj . In this limit P ~ 1 » A;^ f), i.e. r; < 1 

and last term in Eq. (16) can be neglected. As noted above, 
presence of rotational shear in the disc generates SB^f, from 
5Br in this case. 

B. Cyclotron limit: Diffusion in the disc is comparable 
to fluid advection: ^Sv ^ SB^ . In this case fc^ 77 ~ fc^ 2> 1. 
This is the low frequency limit since ct ~ 1 and first 
term in Eq. (16), can be neglected. The low frequency, 
short wavelength dressed ion-cyclotron wave with frequency 
ujc ~ i-L iOH ^ is the normal mode of the system (Pandey & 
Wardle 2008). Therefore, when wave is propagating along 
the vertical field {^J. = 1), the cut-off frequency uic is solely 
determined by the fractional ionization •jjul'^d ~ Ae. In 
PPDs where fractional ionization is small (Wardle & Ng 
1999; Wardle 2007; Bai 2011), left circularly polarized waves 
will have extremely low cut-off frequency. However, the role 
of fractional ionization in determining the cut-off becomes 
redundant if wave in the medium is propagating almost ra- 
dially (/i — >■ 0). In this case, waves will have very low cut-off 
frequency. Therefore, not only the fractional ionization, but 
also the direction of the wavevector determines the cut-off 
threshold of left circularly polarised ion-cyclotron waves. 

The MRI in the cyclotron limit also requires like stan- 
dard MRI picture (a) the generation of radial field due to 
inward fluid motion and (b) the generation of &B^ from 
5Br. However, unlike standard MRI, the azimuthal field is 
not generated by the differential shear of the radial field but 
by the Hall, or, anisotropic (refiected through topological 
switch g) ambipolar diffusion. Clearly, step (b) of the stan- 
dard MRI is now replaced by a totally new step available 
only to diffusive discs. To distinguish this unique feature of 
MRI in diffusive discs, we shall call this case diffusive MRI 
or, DMRI. 

Taking k — > 00, and thus setting a{a) = in Eq. (22) 
or, neglecting the left hand side of Eq. (16), we get following 



^ Tlie term ideal MHD is misnomer here since disc matter is 
partially or, weakly ionised. 

^ The cut-off frequency ldq = '^H^ Eq- (41) given by Kunz (2008) 
correspond to ^ = 1 only. 
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Figure 6. The Diffusive MRI (cyclotron limit, regions II and III) 
and Diffusive Instability (highly diffusive limit, region I, rm > 0) 
of Fig. (2) are highlighted in the above figure. 



dispersion relation 

{jik + flT? Vt) + {fjp -f r?T) 5- + 2 (2 - g) [fjp fjT 



+ r7H + 



= q [ 9 VA + Va + 



(61) 



Since coefficients on the left hand side of preceding equation 
is positive, a positive a requires that the right hand side 
be positive. Thus we arrive at the stability criterion which 
is identical to Eq. (26). This is not surprising since DMRI 
is a subset of MRI and corresponds to region II and III in 
Fig. (6). We have used here Hall-Ohm diffusion dominated 
disc to highlight cyclotron and highly diffusive limits. How- 
ever, Hall-ambipolar dominated disc, e.g. Fig. (4) can as well 
be used for the same pupose. 

The regions II and III in Fig. (6) represent cyclotron 
limit can be directly inferred from marginal stability condi- 
tion [Co ~ 0] also. Thus from Eq. (19) along the Hall axis 
we get 



2 fc2 



1 



m--:, . (62) 

2 - g 

For a Keplerian disc {q = 3/2) in fc — >■ cxj limit, above 
equations describes the interval between rju ~ [—2,-1/2] 
in Figs. (2) and (4). 

Assuming k^ < 1 and thus neglecting k^ terms, for a 
Hall dominated Keplerian disc, from Eq. (61) we get 







1/2 




a ~ 


^ S QuJH 




^k^ilujc 











1/2 



(63) 



Thus the growth rate of ion-cyclotron wave is about one 
and half times the geometric mean of the rotation and ion- 
cyclotron frequencies and attains maximum value only when 
fcz = 1. 

In a purely ambipolar diffusion dominated disc, i.e. set- 
ting rjo = rjH ~ 0, and assuming /x — >■ 0, i.e. when wavenum- 
bers are almost radial, we get 

a^l-gQ, (64) 
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which is identical to the growth rate, Eq. (??). Hence, both 
DMRI and MRI grow at similar rate in ambipolar diffusion 
dominated disc. 

C. Highly diffusive limit: Magnetic Diffusion in the 
disc overwhelms fluid advection; ^5v <^ Jsj . In this case 

Pr? ~ 1 > P, i.e. 77 > 1 and P < 1. The evolution of 
magnetic field which is kinematic in nature in this case is 
solely due etermined by the diffusion. Only first and last 
term in Eq. (16) arc retained. 

In the presence of rotational shear, diffusion of magnetic 
field alone is sufficient to cause the instability in the disc 
(Rudiger& Kitchatinov 2005). It should be pointed out that 
rotation is not essential for this instability. The presence of 
plane shear flow can as well cause this instability. As has 
been noted by Kunz (2008), this is a new kind of instability 
- diffusion instability (DI). 

Unlike MRI or, DMRI, the angular momentum trans- 
port in the disc by DI is not efficient. The physical mecha- 
nism of DI is quite simple: rotational shear generates SB^ 
from 5Br and diffusion completes the loop by generating 
5Br from 6Bj,, i.e. 



SBr SBi 



5B^ 



SBr 



(65) 



It is important to note that although both Hall and ambipo- 
lar diffusion generate SBr from SB^ the generation mech- 
anism in two cases are quite different. Whereas, Hall ef- 
fect through dissipationless diffusion generates SBr from 
5B^, the ambipolar diffusion does it via anisotropic dissipa^ 
tion (Desch 2004). Therefore, suitable disc field and normal 
wavcvector is necessary condition for ambipolar diffusion to 
act like Hall. For purely vertical field, when g = 0, ambipo- 
lar diffusion is incapable of generating SBr from SB^. In this 
case, like Ohm it can only cause the damping of waves. 

The dispersion relation for a Keplerian disc in the highly 
diffusive limit becomes 



&^ + P [2fjo + (1 + M^) tja] <7 + -D = , 



where 



■D = P [??o + (1 + M^) flA + (l- kl 6^^ hi ?7i] 



^{grjA + stth) + [g rjA-s t/h) 



(66) 



(67) 



In the absence of fin, above dispersion relation reduces to 
Eq. (39) of Desch (2004). 

Eq. (66) can be recast in the following form 

{f]p 17t + i?!) P + (''?P + ''7t) + 



= 2^ (gVA + sfju) 



(68) 



from where it is easy to infer general stability criterion for 
a diffusive disc 



gfjA + sfiH> 0- 



(69) 



Clearly excessively diffusive Keplerian disc is unstable if 
the linear combination of ambipolar and Hall diffusion with 
suitable topological coefficients is positive. In the Hall-Ohm 
dominated regime if srjH > the disc becomes diffusively 
unstable. A comparison of Eq. (69) with Eq. (26) shows that 
if the disc is unstable due to DI then it is also unstable due 
to DMRI. However, reverse is not true. This can also be 



seen from Figs. (6). Since DI belongs to region I [beyond 
riH = along the Hall axis, in Fig. (6) corresponding to fi- 
nite ko ,kc, Fig. (1)] guaranteeing that Eq. (26) is satisfied 
in the Hall-Ohm dominated discs. 

In the absence of ambipolar diffusion, as has been shown 
in section 2.2, the dispersion relation, Eq. (66) for Hall-Ohm 
dominated disc can be mapped to a purely vertical field 
threaded disc which has been analysed in considerable detail 
by WSll. Therefore, we shall focus only on Hall-ambipolar 
diffusion dominated disc here. 

The instability in Hall-ambipolar diffusion dominated 
disc will grow at maxirrmm rate 



(To 



{gfjA + sfjn) 



2 [(1-^2)^,2-4^2^)1] 



Va+Vh 



(70) 



which for purely Hall case {rjA = 0) reduces to Eq. (36). 
For purely ambipolar (tjh = 0) case, preceding equation 
becomes 



3 

^0 = - 



2(i + ImI)' 



(71) 



which is similar to Eq. (47). For planar shear gradient 
(dv/dx = 3f2/2) above equation becomes identical to 
Eq. (33) of Kunz (2008). Note that the maximum growth 
rate, Eq. (71) is proportional to the shear gradient and is 
independent of rjA as well as the strength of the background 
magnetic field. However, signature of the magnetic field is 
manifested through the geometric factors g and /u. 

The corresponding most unstable wavevector is also a 
function of g since 



fco VA 



2 Id {1 + Id) fjA 



(72) 



Apart from the shear gradient and geometric factor, the 
wavevector ko also depends on the strength of the ambient 
magnetic field as fco oc 1/B. Thus in the presence of weak 
magnetic field extremely short wavelength fluctuations will 
grow at maximum rate in the ambipolar diffusion dominated 
discs. 

Notwithstanding the independence of maximum growth 
rate on diffusivities in a purely Hall or, purely ambipolar 
case, the general expression Eq. (70), does depend on diffu- 
sivities. In fact, the maximum growth rate, decreases with 
increasing riA/riH- 

In order to see this dependence on the field geome- 
try and wavevector, we plot in Fig. [7(a)-(b)] maximum 
growth rate for two different cases: (a) vary bz and keep kz 
fixed and (b) vary kz and keep hz fixed. We have assumed 
riAlriH = 0.1. The first panel Fig. [7(a)] for fixed kz = 0.9 
suggests that when field is almost vertical, the diffusive in- 
stability grows almost at maximum rate. However, only long 
wavelength fiuctuations, i.e. 2k^riH/3 < 1 can be excited in 
this case. On the other hand in azimuthal field dominated 
disc, fluctuations of smaller wavelength also becomes unsta- 
ble, albeit with reduced growth rate. We see from Fig. [7(b)] 
that the instability will be completely quenched for small kz ■ 
Therefore, in a predominantly azimuthal field threaded disc, 
almost vertical ax;isymmetric fluctuations are most prone to 
diffusive instability. 
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Figure 7. Maximum growth rate versus wavevector in the highly 
diffusive hmit with fixed = 0.9 and varying bz is shown in 4(a) 
and fixed = 0.9 and varying is shown in 4(b). To keep g 
positive, we have assumed = —^/1 — b'j. 



The wave number between [0, kc] where 



kc VA 



3 

2^2 



jflA + Sflu 

Va + Vu 



(73) 



are unstable in the highly diffusive disc. Note that when 
rjH = 0, the cutoff wavelength in AD dominated disc shifts 
to infinity when <? < 0. As a result, AD damps fluctuations 
at all wavelengths. However, when g > 0, only fluctuations 
between [0, kc] can grow, where 



kc va 



(74) 



Choosing maximal value for g — 1/4, i.e. hz = — fe^ = l/\/2, 
and kz = kr = l/\/2, i.e. /.t^/gr = 1, we see from the above 
equation that waves of wavelengths [0, 5 VvaJ will be unsta- 
ble in the AD dominated discs. 

To conclude, excessively diffusive discs can become un- 
stable solely due to interplay between rotational shear and 
magnetic diffusion without significant transport of angular 
momentum. However, the onset of diffusive instability will 
always trigger MRI and thus transport of angular momen- 
tum will always occur in such discs. In discs threaded only 
by a vertical field both DI and DMRI grow at maximum rate 
when the wavevector has only vertical component. However, 
critical wavenumber window increases with decreasing ver- 
tical field. 



5 DISCUSSION 

Analysis of the MRI in three dimensional diffusion space 
for sufficiently general geometry is difficult. Often, proto- 
planetary discs (PPDs) have layered structure with surface 
layer well coupled to the magnetic field whereas midplane 
is magnetically dead (Gammie 1996). Further, Hall-Ohm 
is the dominant diffusion at 1 — 2 scale heights whereas 
ambipolar-Hall dominates at 3 — 4 scale heights from the disc 



midplane (Wardle 2007; Bai 2011). This provides a physi- 
cal motivation to analyse a simplified dispersion relation in 
two dimensional Hall-Ohm or Hall-ambipolar diffusion plane 
which also lends to an easy geometric interpretation. Such 
a representation of the dispersion relation allows us to geo- 
metrically infer the maximum growth rate of the instability. 
However, the stability criterion Eq. (28) pertains to general 
diffusive disc and allows us to infer following conclusion. 



5.1 The diffusive discs are likely to be unstable to 
almost radial fluctuations. 

This can be seen from Eq. (28) which can be written as 
- {m + 2hzv\/Q) 



(^—kr b^^ 



(75) 



Assuming positive left hand side in the preceding equation 
above criterion is easily satisfied for non-zero tja and kz — > 0. 
Therefore, in the presence of almost radial fluctuations, the 
diffusive disc is always unstable provided modes with short 
enough wavelength fits in the disk, i.e. kH >1. 

Note that magnetic tension suppresses MRI growth 
for wavelengths k~^ < va/^ = H/^/]3 (where plasma 
/3 = c^/v^ is the ratio of magnetic to gas pressure) . Only 
when the disc is subthermal (/? > 1), MRI unstable fiuctu- 
ations will fit in the disc. For example, in PPDs, it is quite 
possible that MRI active layer may reside in the upper disc 
when plasma /3 is typically much larger than 1 (Bai 2011). 
Therefore if the upper layer in the inner PPDs (< 10 AU) 
is subthermal, it may become magnetorotationally unstable 
to almost radial fluctuations. At outer disc, even the disc 
midplane may become unstable to such fluctuations. How- 
ever, non-subthermal diffusive discs, MRI modes can not fit 
within the disc scale heights. The disc wind may efficiently 
remove angular momentum in such diffusive discs (Wardle 
& Konigl 1993). 



5.2 Upper and middle layer of PPD is most 
susceptible to MRI and DMRI. 

As has been noted above, owing to its layered structure, 
diffusivities in PPDs vary with height and thus the criti- 
cal wavelength which is a function of diffusivity also varies 
with height. In the upper layer of the disc when Hall- 
ambipolar is the dominant diffusion, fiuctuations of arbitrar- 
ily small wavelength becomes magnetorotationally unstable 
when kz — >■ 0. This can be seen from Eq. (60) which in small 
kz limit gives 



kc VA 



_i_ 



(76) 



Therefore, fcz fcc — > when kz — 0. Since fiuctuations of ar- 
bitrary small wavelength can easily fit within the disc scale 
height (for right range of plasma beta), upper layer of the 
disc belonging to regions II and HI in Hall-ambipolar diffu- 
sion plane [Fig. (4)] will be susceptible to MRI and DMRI in 
the presence of favourable fleld topology, i.e. when B Q, < 
and g > 0. 
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5.3 MRI operates much closer to the midplane of 
PPDs when field is non— vertical. 

In the middle and midplane region of PPDs where Hall- 
Ohm diffusion is the dominant diffusion (Wardle 2007; Bai 
2011), the critical wavenumber is given by Eq. (40). From 
this expression, we infer that in purely Hall dominated discs, 
waves propagating along the magnetic field, (kz = 1) have 
the shortest wavelength. Indeed, setting 770 = in Eq. (40), 
we see that fcz fee — > when — > 0. Therefore, the shortest 
wavelength that fits within scale height is those that are 
propagating vertically. 

In the absence of Ohm diffusion, for fc^ = 1, the criti- 
cal wavenumber Eq. (40) has following dependence on Hall 
diffusivity when field is almost vertical {bz 1) 



fee VA 



1 



(77) 



Q J 2 (j?H + 1/2) " 

We have shown in section 3.1 that in Hall-Ohm dif- 
fusion dominated discs general field topology with oblique 
wavevector can be mapped to pure vertical field and pure 
vertical wavenumber. What is the role of field topology and 
obliqueness of wavevector in such a disc? In order to an- 
swer this question we note from Eq. (40) that when disc is 
threaded prcdoniiuantly by azimutlial field i.e. bz — > crit- 
ical wavenumber can become arbitrarily large since fe^ oc 
1/bz- However, no matter how small Ohm diffusion is, it 
can not be ignored when fi^ qjj term in the denominator of 
Eq. (40) becomes comparable to rio, i.e. bz ~ rjo/VH ~ l//3e. 
Thus for weak vertical field, critical wavelength becomes 

3 bzfju ^31 
2 /j2 j)2 ~ 4 ?7o ■ 



fee VA 

n 



(78) 



How small this vertical field can become? To answer this, we 
recall that the stability condition in a Hall-Ohm dominated 
disc, Eq. (29), for a dustlcss disc can be written as 

Bz > -5.2 10-'''X-^^ f^iyr G . (79) 

where Xe-s = Xe/W~^ and fiiyr ~ ^/ lyr. Thus field can 
not become smaller than m G or else Hall will cease to affect 
MRI. Clearly, coupling of the field to the disc matter sets 
the limit on the smallest vertical field. If the field is smaller 
than Eq. (79), the field is not capable of exciting MRI in the 
disc. 

It is clear from Eqs. (77)-(78) that in Hall- Ohm dif- 
fusion dominated disc, when vertical field is weak, critical 
wavenumber scales with Ohm diffusion as k^ oc l/rjo- Re- 
call that in a purely vertical field threaded disc fe^ oc l/rjo 
(Wardle & Salmeron 2011). This implies that MRI can op- 
erate for much larger values of Ohmic resistivity than is 
possible in a purely vertical field threaded disc. Therefore, 
we may infer that general field topology permits MRI to op- 
erate much closer to the midplane than is otherwise possible 
in discs threaded purely by vertical magnetic field. 

Although our investigation is limited to the linear anal- 
ysis, it is important to connect it to nonlinear investigation 
of the non-ideal regimes in PPDs. The result of numeri- 
cal simulation of PPDs in Hall-Ohm regime (Sano & Stone 
2002) suggests that the quenching of turbulence continues to 
be set by Ohm diffusion even in the presence of Hall effect. 
However, as has been shown by WSll, Sano & Stone (2002) 
assume that Ohm diffusion dominates Hall in their Simula^ 
tion. Therefore, it is not surprising that the results of Sano 



& Stone (2002) shows that MRI turbulence is quenched by 
Ohmic diffusion. The effect of Hall diffusion on the satura- 
tion of the instability is yet to be investigated. 



5.4 Purely vertical field and transverse 
fiuctuation is fundamental to MRI in 
Hall-Ohm dominated disc. 

The general field topology and obliqueness of wavevector 
merely rescales the growth rate and wavevector window. 
However, due to such scaling, fluctuations with smaller 
wavenumber becomes susceptible to MRI in azimuthal field 
dominated discs. Although, discs can become MRI suscepti- 
ble to broader fluctuation spectrum for increasingly weaker 
vertical field, no matter how small the vertical field is, its 
presence is necessary to excite the MRI. Discs with purely 
toroidal field are magnetorotationally stable to axisymmet- 
ric fluctuations. 



5.5 PPDs are also susceptible to purely diffusive 
instability. 

PPDs are not only prone to MRI or, DMRI but can also 
be subject to shear driven diffusive instability. This is due 
to the fact that the magnetic field drift in a weakly ionised 
matter can occur not only due to advective fluid motion but 
can also be caused by the plasma-neutral collision. As a re- 
sult, in the high frequency limit, when field drift is solely due 
to the plasma-neutral collision, waves in a differentially ro- 
tating disc can become unstable due to slippage of the field 
against the sea of neutrals- the diffusion instability (DI). 
On the other hand, in the low frequency limit, when fluid 
advection is completely offset by the field diffusion, disc is 
prone to diffusive MRI (DMRI). The maocimum growth rate 
of both DI and DMRI arc similar. Since DI is high frequency 
counterpart of DMRI it implies that much before the onset 
of DMRI, the DI will destabilise the disc without significant 
transport of angular momentum. Since only weak turbulence 
is excited in excessively diffusive discs [Fig. 13, Bai & Stone 
(2011)], DI may not impact the angular momentum trans- 
port due to DMRI. 



6 SUMMARY 

The MRI in a partially ionised diffusive disc has been anal- 
ysed in the present work. Following is the brief summary of 
the work. 

1. The general stability criterion of a magnetised, dif- 
ferentially rotating diffusive disc [with a rotation profile 
(X r~'' , < q < 2] can be stated purely in terms of 
magnetic diffusivities as 



g'nA+ SIJH > 



2^ v\/Q, 



(80) 



This criterion suggests that subthormal diffusive discs are 
always unstable to almost radial fluctuations. 

2. The MRI can operate much closer to the midplane in 
discs threaded by non-vertical magnetic field. 

3. Whereas discs which are Hall-Ohm stable may or 
may not be ambipolar unstable, Hall-Ohm unstable discs 
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in the presence of non-vertical field and oblique wavevector 
can always become ambipolar unstable. 

4. Favourable field topology and oblique wavevector not 
only make disc magnetorotationally unstable in ambipolar 
regime, but simultaneously enlarges the domain over which 
Hall diffusion can assist MRI. 

5. The vertical magnetic field and transverse axisym- 
mctric fluctuations (vertical wavevector) plays fundamen- 
tal role in destabilising Hall-Ohm diffusion dominated disc. 
Furthermore, MRI grows at maximum rates in vertical field 
threaded discs. The non-vertical field and oblique wavevec- 
tor can be mapped to a purely vertical field and vertical 
wavevector. 

6. The maximum MRI growth rate of the ambipolar dif- 
fusion assisted MRI is two fifth of ideal MHD case (3/4 0) 
for long wavelength fluctuations. Short wavelength fluctusr 
tions are likely to be damped by the ambipolar diffusion. 

7. Excessively diffusive discs can become un.stable solely 
due to interplay between rotational shear of the disc and 
field diffusion. The ensuing diffusive instability grows at the 
same maximum rate as MRI. 
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